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ABSTRACT 

The Lya emission has been observed from galaxies over a redshift span z ~ — 8.6. However, the 
evolution of high-redshift Lya emitters (LAEs), and the link between these populations and local 
galaxies, remain poorly understood. Here, we investigate the Lya properties of progenitors of a local 
L* galaxy by combining cosmological hydrodynamic simulations with three-dimensional radiative 
transfer calculations using the new ART^ code. We find that the main progenitor (the most massive 
one) of a Milky Way-like galaxy has a number of Lya properties close to those of observed LAEs at 
z ~ 2 — 6, but most of the fainter ones appear to fall below the detection limits of current surveys. The 
Lyo; photon escape fraction depends sensitively on a number of physical properties of the galaxy, such 
as mass, star formation rate, and metallicity, as well as galaxy morphology and orientation. Moreover, 
we find that high-redshift LAEs show blue-shifted Lya line profiles characteristic of gas inflow, and 
that the Lya emission by excitation cooling increases with redshift, and becomes dominant at z > 6. 
Our results suggest that some observed LAEs at z ^ 2 — 6 with luminosity of LLya ~ 10^2-43 g-i 
may be similar to the main progenitor of the Milky Way at high redshift, and that they may evolve 
into present-day L* galaxies. 

Subject headings: radiative transfer - line: profiles - hydrodynamics - cosmology: computation - 
galaxies: evolution - galaxies: formation - galaxies: high-redshift 



1. INTRODUCTION 

The Lya emission from young ga laxies can be a pow- 
erful probe of the early universe (jPartridge fc PeeblesI 
^967[ ICharloTXF all 1993). Recent narrow-band deep 
imaging surveys using large-aperture telescopes have 
detected a large number of Lya emitting galaxies, 
or Lyg emitters (LAEs), at redshifts z > 3 (e.g., 
Hu fc McM ahon' ^1996^ 'Cowi e Hul IT998I: iSteidel et al.' 

jMalhotra & Rhoads I2QQ4I: pTa niguchi et al.^ 
2QQ5I : iKashikaw a et al.l 12 006: Shimasak u et al.l 12006 : 



Ive et al.L2q06: Hu fc Cowie 2006: Gronwall et all 

Ouchi et all 120081: iHu et al.l 120101 : lOuchi et al.l 

Lehnert et al.l 120101 ). By combining Lya emission with 
broad-band continuum, multi-wavelength observations 
are beginning to address the physical properties of 
these high-redshift LAEs (e .g., Gawiser et al. 2006; 
Gronwall et al. 2^2^^^^ 2007; Nilss on et al.ll2007 : 
Pirzkal et aL .20071: iLai et al.l i200 8: Ouchi et al." "2008; 
Pentericc i^t all 120091: lOno et al.H 2010a.b: Haves et al 
201Q: iFin kelstein et al.l 120111 : Eflsson & M0ller 2011 
Acquaviv a et al. 2011). It has been suggested that these 
objects are mostly compa ct, young ga laxies with low 
metallicity. In addition, I Ouchi et al.l (j2Q08) studied 
the evolution of equivalent widths (EWs) and the 
characteristic ^Lya with redshift from z ~ 3 to z ~ 6, 
and found that the mean EW increased with redshift, 
while th e did not change significantly. More 

recently, iCiardullo et al.l (|2011f ) studied the evolution 
of luminosity function (LP) and EW from z = 2.1 to 
3.1, and found that increases from 10^^'^ erg s~^ at 
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2007; 



2010; 



z = 2.1 to 10"^^ erg s"^ at z = 3.1. iBlanc et al.l (l2011f ) 
studied the Lya properties of LAEs in the redshift 
range z = 1.9 — 3.8 from the Hobby Eberly Telescope 
Dark Energy Experiment (HETDEX) Pilot Survey, and 
showed that the median Lya escape fraction (/esc) was 
^ 29%, and it does n ot evolve sig nificantly with redshift. 
On the other hand, 'Have s et"al.l (j2011l) suggested that 
/esc monotonically increases between redshift and 6, 
which implies that high-z galaxies tend to be LAEs. 

While high-redshift LAEs have been studied with large 
samples in the redshift range of z ~ 2.2 — 6.6, there is 
only a limited number of observations on LAEs at 2; < 
1. Some local star- forming galaxies have been studied 
by various wavelengths and show a complex structure of 
Lyg an d UV conti n uum (jHaves et al.ll2007l : lOstlin et al.l 
I2009f ). lAtek et all (l2009f ) showed that the /esc of local 
LAEs have a large dispe rsion, ranging from ~ 3 to 100 %. 
In addition, Deharveng et al.l (j^08) studied a sample of 
96 local LAEs at z = 0.2 — 0.35 from UV space telescope 
GALEX^ and found that these LAEs have similar EW 
distrib ution as those at z = 3.1. Recently, iCowie et al.l 
(|2010f ) have studied z ^ 0.3 LAEs with a larger sample, 
and showed that these LAEs are more compact, and have 
lower metallicity than UV-conti nuum selected galaxies at 
the same redshift. In addition, IFinkelstein et al.1 ()2009l ) 
suggested, from fitting of spectral energy distributions 
(SEDs), that low-z LAEs are significantly more massive 
and older galaxies than their high-z counterparts. 

One of the important issues in galaxy evolution is how 
high-re dshift LAEs evolve in to galaxies in the local uni- 
verse. iGawiser et al.l (j2007l ) suggested, from clustering 
analysis, that most z = 3.1 LAEs evolve to present- 
day galaxies of < 2.5 L*, unlike other populations which 



typically evo l ve int o more massive galaxies. Moreover, 
iGuaita et alJ (|2Q1Q| ) indicated that LAEs at z = 2.1 were 
building blocks of present-day L* galaxies such as the 
Milky Way (MW). 

However, the link between high-redshift LAEs and lo- 
cal galaxies, and the probability of these LAEs evolv- 
ing into present-day L* galaxies are not well constrained 
from observations. In order to address these questions, 
one may use the Milky Way as a local laboratory. More- 
over, since Lya emission has been detected from the 
most distant galaxies, understanding of the Lya prop- 
erties of the Milky Way progenitors will provide an im- 
portant clue to the formation of early galaxies. To date, 
there are only a limited numbe r of theoretical stud- 
ies on thi s important issue fe.g.. 'Salvado ri et al.l I2Q1QI : 
iDaval fc L ibeskind 201 ^ Bo th tSalvadori Tt al.l l\201(t ) 
and \D^\ Libeskir^ (|20Tll ) focused on MW progen- 
itors at z ^ 6 constructed from semi-analytical merger 
trees and a cosmological smoothed particle hydrodynam- 
ics (SPH) simulation, respectively. They both used the 
same analytical prescription of Lya emission in which 
the intrinsic Ly g luminosity scales linearly with the star 
formation rate (Daval et al.l I2QQ8] ). However, because 
the Lya properties depend sensitively on a number of 
factors, such as the scattering and propagation of the 
photons in the inhomogeneous medium, the dust content 
of the gas, the ionization structure, the UV continuum, 
and the photon escape fraction. Such a complicated pro- 
cess can only be probed by comprehensive Lya radiative 
transfer calculations combined with realistic simulation 
of galaxy formation. As we will show in this work, our 
detailed Lya modeling on a high-resolution cosmologi- 
cal simulation produce a number of Lya properties such 
as the luminosity functions at different redshifts in good 
agreement with observations. Moreover, in order to in- 
vestigate the evolution of LAEs, we need to study the 
progenitors of the MW at different redshifts systemati- 
cally, not just at a specific time. 

In this paper, we investigate the Lya properties of MW 
progenitors over a wide redshift range of z ~ — 10, 
by combing cosmological SPH simulation of a MW-like 
galaxy from Zhu et al. (in preparation) with 3D RT 
calculations using t he newly improved ART^ code by 
lYajima et al.l (|2Qlll ). Our hydrodynamic simulation in- 
cludes important physics of both dark and baryonic mat- 
ter, and has high resolutions to track the formation his- 
tory of the MW. Our RT calculations include both Lya 
resonant scattering and continuum emission, and are 
done on an adaptive-mesh refinement grid, which cov- 
ers a large dynamical range and resolves the small-scale 
structures in high-density region. Interstellar dust is also 
taken into account to accurately estimate the /esc of Lya 
photons and UV continuum, and the EWs. 

The paper is organized as follows. We describe our 
cosmological simulation in §2, and the RT calculations 
in §3. In §4, we present results of the Lya properties of 
MW progenitors from redshift 10 to 0, which include the 
Lya surface brightness, Lya luminosity, /esc, EW and 
line profile. In §5, we discuss the dependence of /esc on 
physical properties, LAE fraction in our galaxy sample, 
Lya escaping angle and the contribution from excitation 
cooling to Lya emissivity, and we summarize in §6. 

2. GALAXY MODEL 



The cosmological simulation presented here follows the 
formation and evolution of a Milky Way-size galaxy and 
its substructures, as described in detail in Zhu et al. 
(in preparation). The simulation includes dark mat- 
ter, gas dynamics, star formation, black hole growth, 
and feedback processes. The init ial condition is orig i- 
nally from the Aquarius Project (|Springel et al.ll2008f ). 
which produced the largest ever particle simulation of a 
Milky Way-sized dark matter halo. The hydrodynam- 
ical initial condition is reconstructed from the original 
collisionless one by splitting each original particle into 
a dark matter and g a s par ticle pair, as adopted from 
IWadepuhl fc Springell (|2QTTh . 

The whole simulation falls in a periodic box of 
100 h~^Mpc on each side with a zoom-in region of a size 
5x5x5 h~^Mpc^. The spatial resolution is ~ 250 
pc in the zoom-in region. The mass resolution of this 
zoom- in region is 1.8 x 10^ h~^MQ for dark matter par- 
ticles, 3 X 10^ H-^Mq for gas, and 1.5 x 10^ h-^M^ 
for star particles. The cosmological parameters used in 
the simulation are Vtrn = 0.25, I^a = 0.75, erg = 0.9 
and h = 0.73, consis tent with the five-year results of the 
WMAP (|Komatsu et a l. 2Q09|). The simulation evolves 
from z = l21 to z = 0. 

The simulation was performed using the parallel, N- 
body/SPH code GADG ET-3, which is an imp r oved ver- 
sion of that described in iSpringel et al.l ()2001f ): iSpringell 
(2005.). For the computation of grav itational f orces, 
the code uses the "TreePM" method (IXul 112951) that 
combines a "tree" algorithm (jBarnes &: HutI [19861 ) for 
short-ran ge forces and a Fourier t ransfo rm particle- mesh 
method (jHocknev fc E astwood! Il981f ) for long-range 
forces. GADGET iniplement s the entropy -conserving 
formulation of SPH (|Springel Hernauistl I2002D with 
adapt ive particle smoothing, as in iHernquist fc Katzl 
(|1989f ). Radiative cooling and heating processes are 
calculated assuni i ng collision al ionization equilibrium 
(jKatz et all 119961 : iDave et al.l [r999). Star formation is 
modeled in a multi-phase ISM, with a rate that follows 
the S chmidt-Kennicutt Law (Schmidt 1959; Ke nnicutFI 
I1998D . Feedback from supernovae is captured through a 
multi-phase model of the IS M by an effective equation o f 
state for star- forming gas (|Springel & HernQuistll2003l) . 
The UV background model of iHaardt fc Madau I (|1996f ) 
is used. 

Black hole growth and feedback are al so included in 
our simulation based on the m odel of ISpringel et al.l 
(j2005b[ ): iDi Matteo et al.l (j2005f ). where black holes are 
represented by collisionless "sink" particles that inter- 
act gravitationally with other components and accrete 
gas from their surroundings. The accretion rate is esti- 
mated from the local gas d ensity and sound speed us- 
ings; a spherical Bondi-Hoyle (|B"ondilll952l : lBondi fc Hovld 
119441 : iHovle fc Lvttletonlll94ir ) model that is limited by 
the Eddington rate. Feedback from black hole accretion 
is modeled as thermal energy, ^ 5% of the radiation, in- 
jected into surroundin g ga s isotropically, as described in 
ISpringel et all(j2005br ) and lDi Matteo et al.l (j2005D . This 
feedback scheme self-regulates the growth of the black 
hole and has been demonstrated to successfully repro- 
duce many observed prope r ties of local elliptica l galaxies 
fe.g.. ISpringel et al.ll2005al : iHopkins et al.ll2006D and the 
most distant quasars at z ^ 6 (|Li et al.l 120071 ). We fol- 
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low the black ho l e seed ing scheme of iLi et aTl (|2QQ7[ ) and 
iDi Matteo et alJ (|2QQ8l ) in the simulation: a seed black 
hole of mass Mbh = 10^ H'^Mq was planted in the grav- 
itational potential minimum of each new halo identified 
by the friends-of-friends (FOF) group finding algorithms 
with a total mass greater than 10^^ h~^MQ. 

The galaxies in each snapshot above redshift are 
building blocks of the MW, a present-day L* galaxy. We 
therefore define them as "progenitors" of the MW, and 
the most massive progenitor at any given timestep as 
the "main progenitor" . In this paper, we explore the Lya 
properties of ^ 60 most massive progenitors of each snap- 
shot for 15 snapshots in the redshift range z = — 10.2, 
which gives a total sample of 941 galaxies. 

3. RADIATIVE TRANSFER 

The RT calculations are done using the 3D Monte 
Carlo RT code. All- wavelength Radiative Transfer with 
Adaptive Refinement Tree (ART^), as recently developed 
by lYaiima et all (1201 ID. AR T^ was improved over the 
original version ofE i et al.l ()2008l ). and features three 
essential modules: continuum emission from X-ray to 
radio, Lya emission from both recombination and col- 
lisional excitation, and ionization of neutral hydrogen. 
The coupling of these three modules, together with an 
adaptive refinement grid, enables a self-consistent and 
accurate calculation of the Lya properties, which depend 
strongly on the UV continuum, ionization structure, and 
dust content of the object. Moreover, it efficiently pro- 
duces mult i- wavelength properties, such as the spectral 
energy distribution and images, for direct comparison 
with multi-band observations. The detail ed imple menta- 
tion s of the ART^ code are described in iLi et al.l (|2Q08i) 
and lYajima et all ()2011f ). Here we focus on the Lya cal- 
culations and brieffy outline the process. 

The Lya emission is generated by two major mecha- 
nisms: recombination of ionizing photons and collisional 
excitation of hydrogen gas. In the recombination process, 
we consider ionization of neutral hydrogen by ionizing 
radiation from stars, active galactic nucleus (AGN), and 
UV background (UVB), as well as by collisions by high- 
temperature gas. The ionized hydrogen atoms then re- 
combine and create Lya photons via the state transition 
2P IS. The Lya emissivity from the recombination is 

C = faC^Bhi^arienmi, (1) 

where ae is the case B recombination coefficient, and 
fa is the average number of Lya photons produced 
per case B rec ombination. Here we use ae derived in 
iHui & GnedinI ([T997i ). Since the temperature depen- 
dence of fry is not strong, = 0. 68 is assumed every- 
where (jOsterbrock fc Ferlandll2QQ6f ). The product hua is 
the energy of a Lya photon, 10.2 eV. 

In the process of collisional excitation, high tempera- 
ture electrons can excite the quantum state of hydrogen 
gas by the collision. Due to the large Einstein A co- 
efficient, the hydrogen gas can occur de-excitation with 
the Lya emission. The Lya emissivity by the collisional 
excitation is estimated by 

e'f' = CLycnenm, (2) 

where Ci^ya is the collisional excitation coefficient, 
CLya = 3.7 X 10-^^exp(-/iz/^/fcT)T-i/2 g^gg g-l ^^3 



(jOsterbrock Ferlandll2006D . 

Once the ionization structure have been determined, 
we estimate the intrinsic Lya emissivity in each cell by 
the sum of above Lya emissivity, = e^®^ + e^^^. 

In RT calculations, dust extinction from the ISM is in- 
cluded. The dust content is estimated according to the 
gas content and metallicity in each cell, which are taken 
from the hydrodynamic simulation. The dust-to-gas ra- 
tio of the MW is used where the metallicity is of Solar 
abundance, and it is linearly interpolated for other metal- 
licity. We u se the stellar population synthesis model of 
GALAXEV (jBruzual fc Charlotll2003l ) to produce intrin- 
sic SEDs of stars for a grid of metallicity and age, and 
we use a simple, broken p ower law for the AGN ()Li et al.l 
I2OO8I ). A ISalpeterl ([T955f ) initial mass function is used in 
our calculations. 

In this work, we apply ART^ to the 60 most massive 
progenitors of each snapshot for 15 snapshots at redshifts 
spanning z = — 10.2. In our post-processing procedure, 
we first calculate the RT of ionizing photons (A < 912 A) 
and estimate the ionization fraction of the ISM. The re- 
sulting ionization structure is then used to run the Lya 
RT to derive the emissivity, followed by the calculation 
of non-ionizing continuum photons (A > 912 A) in each 
cell. Our fiducial run is done with TVph = 10^ photon 
packets for each ionizing, Lya, and non-ionizing compo- 
nents. Because the spatial resolution of the cosmological 
simulation is not adequate to resolve the multiple phase 
of the ISM, we assume a single-phase medium in each 
density grid. The highest refinement of the grid corre- 
sponds to a cell size comparable to the spatial resolution 
of 250 pc in comoving coordinate of the hydrodynamic 
simulation. 

4. RESULTS 

4.1. The Formation History of The Milky Way 

Figure [T] shows the evolution of the MW galaxy from 
redshift z=10 to z=0 from the cosmological simulation. 
The gas follows the distribution of dark matter in fila- 
mentary structures, and stars form in high density re- 
gions along the filaments. The most massive galaxy re- 
sides in the intersection of the filaments, the highest den- 
sity peak in the simulated volume where gas concentrates 
in the deep potential well. The MW galaxy is formed by 
gas accretion and merger of subhalos, it has the last ma- 
jor merger at redshift z ^ 2 (Zhu et al, in preparation). 

Stars start to form at z ~ 15 by accretion of cold gas. 
As shown in Figure [21 the star formation rate (SFR) 
of the main progenitor reaches ~ 10 Mq yr~^ at z ~ 
10, and it peaks at ~ 58 M© yr~^ at z = 5.2, owing 
to merger of gas-rich protogalaxies. Galaxy interaction 
induces gravitational torques and shocks, which trigger 
global starburst. After that, the SFR generally decreases 
except a boosted bump at 2; ~ 2 when the last major 
merger takes place. 

Shown also in Figure [2] are the SFRs of the top 60 
progenitors at different redshifts. The median and total 
value of these SFRs gradually increase until it reaches the 
peak at z ~ 3.5, after which it decreases rapidly by over 
an order of magnitude. This evolution is in broad agree- 
ment with the observed cosmic star formation history 
(|Hopkins fc Beacomll2006l ). although the simulation box 
is somewhat small to discuss such a statistical property. 
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Fig. 1. — Evolution of the MW galaxy from redshift z=10 to z=0. The images are projected density of gas and stars. The box size is 1 
Mpc in comoving scale. The gas follows the distribution of dark matter and shows filamentary structures. Stars and galaxies form in high 
density regions along the filaments. 
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Fig. 2. — The growth history of the MW galaxy. Top panel: the 
individual star formation history of the most massive 60 progeni- 
tors at each redshift. Red filled circles represent the main progen- 
itor (the most massive one), and the blue filled circles indicate the 
median value of our sample galaxies represented by black open cir- 
cles. Bottom panel: the accumulated mass of the main progenitor 
at different redshift. Filled circle represents the total mass, while 
the open circle represents the stellar mass. 



The star formation at high redshifts [z > 4) is largely fu- 
eled by inflow of cold gas and mergers of gas-rich halos, 
while the rapid decline of SFR at z < 2 is mainly caused 
by feedback from stars and AGNs, and the depletion of 
cold gas. 

The main progenitor has a total mass of ~ 2.6 x 
10^° M0, and a stellar mass of - 8.5 x 10^ at z - 10. 
It evolved into a disk galaxy at 2; = with a total mass of 
1.6 X 10^^ M0 and stellar mass of ~ 10^^ M©, as observed 
in the MW galaxy. 



4.2. Lya Surface Brightness 

Figure [3] shows the Lya surface brightness of the MW 
galaxy at different redshifts, contrasted with distribu- 
tions of gas and stars of the galaxy. To facilitate compar- 
ison with observations, we adopt an intensity threshold. 
Shy a = 10 ergs s~^ cm~^ arcsec"^, fro m a recent sur- 
vey of extended Lya sources at 2; ^ 3 by iMatsuda et al.l 
(|2011f ) to show the contours. The Lya distribution ap- 
pears to trace that of the gas. At z > 7.2, the galaxy 
is small and compact, and the Lya emission is confined 
in the central high-density region. As the galaxy grows 
in mass and size, the Lya emission becomes more ex- 
tended. At 2: = 4.2 — 3.1, the gas structure is irregular 
due to infall along with filament of the main halo. At 
z ~ 0.0, the galaxy shows a disk geometry with spi- 
ral structures. Indeed, the Lya map shows filamentary 
structures at high redshift, and spirals at z ~ 0.0. We 
note that some of the extended Lya so urces in the recent 
observations by IMatsuda et al.l ()2011l ) , which are called 
Lya blobs, show filamentary structures. Our galaxy 
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Fig. 3. — Evolution of the Lya surface brightness of the MW galaxy with redshift. The first and second columns show the distribution 
of gas (in green) and stars (in blue) from the cosmological simulation, respectively. The third column shows the surface brightness of Lyo: 
in log scale with units of ergs s"-*^ cm~^ arcsec"^. The box size is 200 kpc in physical coordinates. Note the galaxy at z = is artificially 
placed at 2; = 0.1 in order to obtain the luminosity distance. The right column shows the contour of Lya surface bri ghtness at th e level 
of lO"-*^^ ergs s"-*^ cm~^ arcsec"^, similar to the threshold of recent observations of extended Lya sources bv Matsud a et al.l (|201in . The 
meshes are coarse grained to be 4 kpc in physical scale to make the smooth contours. 
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Fig. 4. — Radial distribution of the surface brightness of the MW 
galaxy at different redshift. The solid line shows the Lya surface 
brightness in units of lO"-*^^ ergs s"-*^ cm~^ arcsec"^, while the 
dotted line shows the surface brightness of UV continuum at 1450 A 
in rest frame in units of 10~^ ergs s"-*^ cm~^ Hz"-*^ arcsec"^. 
The surface brightness distribution is derived from the projected 
flux images in flgure[3]by taking the mean surface brightness as a 
function of the distance. 

at z = 3.1 has a Lya distribution of ~ 50 kpc with 
a surface brightness above the observational threshold. 
However, the size is smaller than the observed giant 
Lya blobs of > 100 kpc. Such large blobs are prob- 
ably produced by systems of Mhaio ^ 10^^ M©, more 
massive than our model. In addition, it was suggested 
that ex tended Lya sour ces of > 40 kpc become rare at 
z < 1 (jKeel et al.l l2009f ). Recent UV surveys detected 
Lya from local star- forming galaxies at 2; < 1, which 
showed ~ 10 kpc Lya distribution above t he threshold 
of Si^ya ~ 10"^"^ ergs s~^ cm~^ arcsec"^ (jHaves et al.l 
120071 : lOstlin et al'll2009f ). With such a detection sensi- 
tivity, our model galaxies at z < 1 show the size of < 20 
kpc, in broad agreement with observations. 

To examine the difference in distribution between stars 
and Lya emission more quantitatively. Figure H] shows 
the surface brightness of UV continuum, which traces 
young stars, and Lya as a function of the distance from 
galaxy center at different time. At high redshift ^ 5), 
the Lya distributes more extendedly than the stars, as 
the UV continuum decreases steeply around the virial 
radius, but at z < 3, both emissions appear to have sim- 
ilar radial distribution. Such a transition is mainly due 
to the difference in Lya production at different epochs. 
The Lya emission is dominated by collisional excitation, 
which depends strongly on the gas density, at z > 5. At 
a later time, Lya from recombination of ionized gas by 
stellar radiation becomes more important, so Lya emis- 
sion follows that o f stars. 

It was shown bv lCowie et al.l ()2011f ) that nearby LAEs 
(z < 1.0) have a variety of morphologies, some are disky. 



while others are mostly compact galaxies. Galaxy mor- 
phology is closely tied to galaxy formation and evolu- 
tion, and it is related to the galaxy mass. As demon- 
strated in Figure [5l which shows a sample of four galax- 
ies of different masses at z = 0.2, the Lya morphol- 
ogy changes with galaxy mass. In galaxies with lower 
mass (Mtot ^ 10^^ ^0)7 the Lya luminosity is low 
{Lhya ^ 10^^ erg s~^), the Lya morphology is highly 
compact. At higher mass (Mtot ^ 10^^ M©), the Lya 
luminosity is high (i^Lya ^ 10^^ erg s~^), Lya morphol- 
ogy shows disky and spiral structures. This plot suggests 
that the various morphologies observed in low-redshift 
LAEs may reflect a wide range of galaxy mass in the 
sample. 

We note that in Figure IH the azimuthally averaged 
surface brightness at z > 3 is somewhat fainter than 
the detection threshold of recent n arrow-band surveys , 
- 10-^^ ergs s-^ cm-2 arcsec-^ (jOuchi et al.l l2008f ). 
However, since the local Lya distribution is inhomoge- 
neous and anisotropic, as shown in Figure [3l so bright 
regions with flux above the threshold may be detectable 
by such surveys. 

The detectability of these galaxies depends strongly 
on the sensitivity of the surveys. In the present work, 
unless noted otherwise, the Lya luminosity is computed 
by collecting all escaped photons without a flux limit. If 
a detection limit of a given instrument is imposed, the 
luminosity of individual galaxies, in particular that of the 
faint ones, may be reduced significantly, as suggested by 
IZheng et al.l (|20Tol ). 

4.3. Evolution of Spectral Energy Distribution 

The corresponding multi-wavelength SEDs of the 
galaxy sample in Figure [3] are shown in Figure [6l The 
shape of the SED changes significantly from z = 10.2 
to z = 0, as a result of changes in radiation source and 
environment, since the radiation from stars, absorption 
of ionizing photons by gas and dust, and re-emission by 
the dust evolve dynamically with time. The Lya line 
appears to be strong in all cases. The deep decline of 
Lyman continuum (/ < 912 A) at high redshifts {z > 7) 
is caused by strong absorption of ionizing photons by the 
dense gas. Galaxies at lower redshift have a higher floor 
of continuum emission from stars and accreting BHs, a 
higher ionization fraction of the gas, and a higher in- 
frared bump owing to increasing amount of dust and 
absorption. Moreover, due to the effect of negative k- 
correction, the flux at observed frame A > 500 /im stays 
close in different redshifts. 

Our calculations show that the main progenitor has a 
flux of fiy = 0.057 mJy at z ~ 6 and 0.02 mJy at z ~ 8.5 
at 850 /im in observed frame. The new radio telescope, 
Atacama Large Millimeter / suhmillimeter Array (ALMA) 
may be able to detect such galaxies at z ~ 6 with ~ 2 
hours inte gration, and at z ^ 8.5 for ^ 20 hours with 16 
antennas (jYajima et al.ll20l'T| ). 

4.4. The Lya Properties 

The resulting Lya properties of the 60 most massive 
progenitors from selected snapshots, and their evolution 
with redshift are shown in Figure H The top panel shows 
the emergent Lya luminosity ^Lyo;- The main progen- 
itor has a luminosity of I/Lya ~ 10^^ erg s~^ at z < 2, 



6 



^ in % 
od 00 2 



C7> CT> CJ> 
O O O 



00 

00 ^ 

II II 



o es 



C7> tn 



Q> o> o> 
o o o 



o 



CN 



II 1 



q Q 



o 

0> Q> 0> 

o o o 







O 
CN 


d 




gos)=1 


II 






II 


Sto 


_l 














_o 


O 


O 




■■sSS^?^-j-.- 



' ' ' I - ' ■ V- ' - r v- ' ^ ■ ' ■ I 






star 



I ■ ■ ■ I 




I ■ ■ ■ I 



1 



Lya 



■ ■ ■ - J- . 





u 

CD 
0) 
U 



E 
o 

'co 
cn 

(D 
CO 

o 



-16 
-17 

-18 

-19 

-20 
U-21 



Fig. 5. — The relation between Lyck; morphology and galaxy mass, exemplified by a sample of four galaxies at z = 0.2 with different 
masses, Mtot = 4.2 x 10^, 1.7 x 10^°, 1.2 x 10^^ and 1.6 x 10^^ M©, respectively. The left and middle columns show the distribution of 
gas (in green) and stars (in blue) from the cosmological simulation, respectively. The right column shows the Lya surface brightness of the 
corresponding galaxies, with resulting luminosity I/Lya = 2.6 x 10^*-*, 9.8 x 10^*-*, 1.8 x 10^^ and 9.8 x 10^^ erg s"-*^, respectively. The box 
size is 250 kpc in physical scale. Note the Lya luminosity is computed by collecting all escaped photons. 
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Fig. 6. — The SEDs of the most massive galaxy in each snap 
shot. Only the galaxy at z = 0.0 is artificially set at z = 0.1, i.e., 
the fiux is estimated by using the luminosity distance of z = 0.1. 



then increases to Li^ya ~ 10 erg s~ at z = 2 — 6, 
owing to the increase of SFR. At high redshift z > 6^ 
the Li^ya decreases to ~ 10^^ erg s~^ due to low escape 
fraction and absorption by dust. At high redshift, the 
galaxy size is small and most of the stars form around 
the galaxy center. Hence, the dust compactly distributes 
around young stars, which effectively absorbs the ionizing 
photons. As a result, the intrinsic Lya emissivity drops 
even though SFR is enhanced by the accretion of cold 
gas. However, most of other model galaxies show that 
the escape fraction of Lya and UV continuum photons 
monotonically increases with redshift because of lower 
dust content. The escaping process of continuum pho- 
tons will be discussed in detail in a forthcoming paper 
(Yajima et al. in preparation). 

Most of the galaxies at z > 3 in this simulation have 
a luminosity below 10^^ erg s~^, the lower limit of many 
current observations using narrow-band filters, so they 
may not be observable. Howe ver, the impro ved sen- 
sitivity of deep survey of Cassat a et al.l (|2011l ) reaches 
Lhya ^ 10^^ erg s~^atz~ 2 — 6.6, which may detect 
more faint galaxies as the ones in our sample. As men- 
tioned in previous sections, the luminosity is calculated 
as the sum of all escaped photons. If we consider only 
pixels brighter than S = 10~^^ ergs s~^ cm~^ arcsec"^, 
the detection threshold of h i^h-redshift LAE surveys 
fe.g. JOuchi et al.ll2008L I2010D . the total flux of the main 
progenitor is reduced by a factor of a few, resulting in 
Li^ya - 1 - 2 X 10^^ erg s"^ in redshift z - 2 - 6, 
which is close to the observed L^^^of LAEs in this 
redsh i ft span (e.^. , iGawiser et al.l 120071: iGronwall et al.l 
[2QQ7I : lOuchi et al.l 120081 : ICiarduno et al. I I2011D . This 
main progenitor has a halo mass of ^ 10^ ''^ Mq at z ^ 6, 
and ~ 10^^ Mq at z ~ 2, in good agreement with sug- 




z 



Fig. 7. — Evolution of the Lya properties of the 60 most massive 
progenitors from selected redshifts. From top to bottom panel is 
the emergent unfiltered Lyo; luminosity (computed by collecting 
all escaped photons), escape fraction of Lya photons over whole 
solid angle, and equivalent width in rest frame, respectively. The 
red filled circle represents the value of the main progenitor at each 
redshift, while the blue filled circles indicate the median value of 
the sample galaxies. 

gestions from clustering analysis by lOuchi et al.l (|2010f ) . 
These results suggest that some the observed LAEs at 
2: ~ 2 — 6 may be similar to the main progenitors of 
MW-like L* galaxies at high redshifts. 

The calculated escape fractions of Lya photons (/esc) 
are shown in the middle panel of Figure [71 Here, the /esc 
is estimated by correcting all escaped photons over whole 
solid angle and dividing by intrinsically emitted photon 
number. Unlike the SFR, the /esc has higher values at 
lower redshift z <2, then decreases gradually to 20 %. 
At z > 4 the /esc increases again. The /esc of the main 
progenitor fluctuates in the range of ~ 20 — 60 %. The 
median /esc at 2 < z < 4 is ~ 30 %, which is consistent 
with the recent observ ation by the HETDEX pilot sur- 
vey (jBlanc et al.ll201ll ). At lower redshift z < 1, there 
is a large dispersion in /esc, similar to the recent obser- 
vation by lAtek et all (| 20091 ) . This large scattering may 
be caused by variation in a number of physical proper- 
ties such as SFR, metallicity, and disk orientation. We 
will discuss the dependence of /esc on these properties in 
detail in Section ISTTl 

We note that the Lya RT calculations in our work, 
which take into account local ionization structure and 
inhomogeneous density distribution of gas and dust, 
produce a smaller escape fraction (/esc ^ 20 — 80 % 

at z ~ 6) than that in previous semi-analytical 

work of iS alvadori et aL| ([2010h (/esc > 80%) and 
iDaval Libeskind (,20nh (/esc - 60 - 90 %), in which 



8 



a uniform slab model was assumed. We find that more 
than half of Lya photons can be absorbed, because dense 
gas and dust around the star-forming and Lya-emitting 
regions absorb the photons effectively. 

The EW of Lya line is defined by the ratio between 
the Lya flux and the UV flux density /uv in rest frame, 
where the mean flux density of A = 1300 — 1600 A in 
rest frame is used. The resulting Lya EWs are shown 
in the bottom panel of Figure [71 Most of the galax- 
ies have E W > 20 A, the y are therefore classified as 
LAEs fe.g.. iGronwan et al.ll20()7l ). The median EW in- 
creases with redshift, from ^ 30 A at redshift z = to 
^ 820 A at z ^ 8.5. This trend is in broad agreement 
with observations that galaxies at higher redshifts ap- 
pear to have hi gher EW than t heir c ounterparts a t lower 
redshifts fe.g.. iGronwah et al.l [20071 : lOuchi et~alll2008[ ). 
The high EW at 2; > 6 is produced by excitation cool- 
ing, which enhances the Lya emission at high redshift, 
but at low redshift it reduces the EW as th e stellar pop- 
ulation ages (e.g.. [Finkelstein et al.[[2QQ9[ ). Recent ob- 
servations of LAEs at 2; ^ 0.2 — 0.4 shows that most 
local LAEs, unlike those at 2; > 3, have EWs less than 
100 A ([Deharveng et al.[[2008[ : [Cowie et al.[[2011[ ). consis- 
tent with the trend seen in our model. 

We should point out that the results presented in Fig- 
ure [3 are "unfiltered" by detection limit, and that we 
caution against taking these numbers too literally when 
compared with a particular survey, because the observed 
properties depend strongly on the observational thresh- 
old. Note also in the current work, we did not include the 
transmission in intergalactic medium (IGM). The Lya 
properties can be changed by IGM extinction. The neu- 
tral hydrogen in IGM at high redshift can scatter a part 
of Lya p hotons, and d e crease the Li^y^ and EW. For 
example, [Laursen et al.[ (|2011[ ) suggested that the IGM 
transmission could be ^ 20 % at z = 6.5. The trans- 
mission depends sensitively on the viewing angle and the 
environments of a galaxy, as it is affected by the inho- 
mogeneous filamentary structure of IGM. 

4.5. The Lya Line Profiles 

The emergent Lya emission line of the main progenitor 
is shown in Figure [51 The frequency of the intrinsic Lya 
photon is sampled from a Maxwellian distribution with 
the gas temperature at the emission location in the rest 
flame of the gas. Our sample Lya lines show mostly 
single peak, common profi les of LAEs obser ved both at 
high redshift {z ^ 6) (e.g., [Quchi et al.[ 2010[ ) and in the 
nearby universe fe.g.. lCowie et al.[[2010[ ). 

In a static and optically thick medium, the Lya pro- 
file can be double peaked, but when the effective op- 
tical depth is small due to high relative gas speed or 
ionization state, there niight be only a single peak 
([Zheng fc Miralda-Escudl[2002[ ). In our case, the flow 
speed of gas is up to ^ 300 km/s, and the gas is highly 
ionized by stellar and AGN radiation, which result in a 
single peak. 

In the case at high redshift 2; > 6, the gas is highly con- 
centrated around the galaxy center, hence they become 
optically thick and cause the Lya photons to move to the 
wing sides. In addition, the profile at 2: = 10.2 shifts to 
shorter wavelength, and it shows the characteristic shape 
of gas inflow (^Zheng fc Miralda-Escudd.2002i) . Although 
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Fig. 8. — The Lya line profile of the main progenitor at different 
redshifts. The dotted- and solid lines are the intrinsic and emergent 
Lya profile, respectively. 



our simula tion includes feedback of stellar wind similar 
to that of [Springel et al.[ (|2005b[ ). the Lya line profile 
indicates gas inflow in the galaxy. Our result suggests 
that high-redshift star-forming galaxies may be fueled 
by efficient inflow of cold gas from the filaments. We 
will study this phenomenon in detail in Yajima et al. 
(in preparation). On the other hand, it was suggested 
that the asymmetrically shifted profile to the red wing 
in so me LAEs can b e made by outflowing gas distribution 
fe.g.. [Mas-Hesse et a l 2003). The growing hot bubble gas 
around star-forming region from supernovae or radiative 
feedback can cause outflowing, neutral gas-shells, which 
result in re d-shifted line profiles. 

Recently, [Yamada et al.[ (|2012[ ) observed a sample of 
91 LAEs at z = 3.1, about half of which show double 
peaks of strong-blue and weak-red features thought to 
be caused by gas outflow, while others show a symmetric 
single peak in which the flux ratio of blue wing to red one 
is about unity. While our model may explain the latter, 
the missing outflow features in our line sample is proba- 
bly due to the hmitations of our current simulations, e.g., 
insufficient spatial resolution and simplified treatment of 
supernovae feedback. 

In addition, the line profiles of galaxies at high redshift 
z > 6 may be highl y suppressed a nd changed by scat- 



y scat- 
][2003; 



tering in IGM (e.g., Santos 2004; ' Diikstra et al. ._ „ 

Zheng et al. 2010; Laursen et al. 20li[), because the Lya 
transmission through IGM is very low at the line cen- 
ter and at shorter w avelengths by the Hubble flow (e.g., 
[Laursen et al.[ [201l[ ). Even at lower redshift z ^ 3, 
the optical depth of IGM can be high depending on 
the viewing angle a nd the location of the galaxy (e.g., 
[Laursen et al.[[2Qll[ ). Therefore, the Lya flux with in- 
flow featured in our model galaxies may be suppressed 
and the shape may change to a single peak with only the 
red wing, or a do uble peak with s trong -red and weak-blue 
as in Figure 7 of [Laursen et al.[ ()2011f ). 
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5. DISCUSSION 

5.1. Dependence ofLya Properties on Galaxy 
Properties 

As shown in previous Sections, the Lya properties vary 
significantly in different galaxies. Here we explore the de- 
pendence on a number of physical properties of a galaxy. 
Figure [9] shows the dependence of escape fraction /esc 
(top panels) and Lya luminosity Li^ya (bottom panels) 
on the galaxy mass, SFR, and metallicity Z. We apply 
a least-absolute-deviations fitting to the data using a 
power-law function, logF = alogX + f3. 

The mass dependence of /esc has a large dispersion, but 
from our fitting, a ^ —0.02 and /3 ^ —0.17, which sug- 
gests that /esc roughly de creases with the total mass, con- 
sistent with the results of lLaursen et al.l (|2009f ). At M ~ 
10^°-^^ Mo, the /esc is mostly constant at - 10 - 30 %. 
In contrast, /esc is more tightly correlated with the SFR, 

with a r 0.08 and /3 r 0.41. At high SFR, dust 

can be enriched quickly by type II super novae, and can 
effectively absorb the Lya photons. In addition, galax- 
ies with high SFR have more hydrogen gas. The gas 
decreases the mean free path of Lya photons, resulting 
in the increase of the dust optical depth which reduces 
the escape fraction. In addition, the /esc decreases with 
metallicity, a ~ —0.35 and /3 ~ —0.65. Since the dust 
content linearly increases with metallicity in our model, 
the Lya photons can be absorbed effectively by gas with 
high metallicity. Th is trend is consistent with observa- 
tional indication by lAtek et aTl (j^09) and Haves et al. 
(|20T0h . 

On the other hand, the luminosity L\^ya has different 
relationships with these properties from the /esc- The 
I/Lya is also roughly correlated with the mass, LLya — 
-j^Q37.7 ^ M^^^ ^ with a large dispersion. Only massive 
galaxies with Mtot ^10^^ have the Lya luminosity 
of LLya > 10^^ erg s~^. This is consistent with sug- 
gestions from clustering analysis o f observed LAEs at 
^ = 2 - 3 (e.g., iGawiser et £01120071 : iGuaita et ani2010[ ). 
In our model, the massive galaxies at 2; = 2 — 3 evolve 
into L* galaxies at z = 0. Hence, o ur results support 
the suggestion bv IGawiser et al.l (|2007[ ) that the observed 
LAEs with LLya > 10"^^ erg s"^ at z = 2 - 3 are likely 
progenitors of local L* galaxies. 

The ^Lyo; has the tightest correlation with SFR among 
the properties investigated here: ^Lya ^ 10^^'^ x SFR^-^^ 
In the literature, a simple linear relation is commonly 
used, with i^Lya (erg s-^) = 1.1 x 10"^^ x SFR (M© yr-^), 
assuming that i^Lya/^Ha = 8.7 (case B). However, our 
result suggests that the relation between LLya and SFR 
becomes somewhat shallower due to the dependence of 
/esc on SFR. Finally, the emergent LLya does not show 
a strong dependence on metallicity, L\^ya — 10^^*^ x 
(Z/Z©)"^-^^. This is due to the fact that, although the 
intrinsic L\^ya increases with halo mass (so does SFR 
and metallicity), the /esc decreases with metallicity, so 
^Lya of higher-metallicity galaxies is suppressed by dust 
absorption. 

We should point out that the large scatter in the cor- 
relations in Figure [9] may be due to the small volume of 
our simulation and the small number of our galaxy sam- 
ple. In addition, as we discuss in Section 15. 6[ a number 
of limitations of our model, such as the simplified ISM 



model and insufficient resolutions, may contribute to un- 
certainty in these relations. Moreover, the luminosity 
scaling relations may change under some specific detec- 
tion limits. We will study these relations in detail with 
improved model and simulations in future work. 



5.2. Redshift Dependence of LAE Fraction 

The number fraction of LAEs (/lae) in our sam- 
ple is shown in figure [TOl The detection limit of Lya 
varies in different surveys. At high redshift s {z > 
3), the LAE detection in most of observations has 
been confined to Li^y^ ^ 10^^ erg s~^. Here, we de- 
rive the /lae with three Li^y^ thresholds, Li^y^ ^ 
10"^^ 10^\ 10^'^ erg s-i with EW of > 20 A. The /lae 
with Li^ya ^ 10^^, 10^^ erg s~^ rapidly increases from 
z = to ~ 5, and then remains nearly constant with 
higher values > 0.8 at z > 5. The trend is roughly sim- 
ilar to the SFR history (Figure [2]). Since the LLya is 
tightly correlated with the SFR (Figure [9]), the number 
of galaxies with Li^ya ^ 10^^, 10^^ increases dit z ^ — 4. 
On the other hand, although the SFR decreases at z > 4, 
the /esc increases due to low metallicity. Hence, the /lae 
does not decrease at 2: > 4. 

Meanwhile, the /lae with Li^y^ ^ 10^^ erg s~^ is 
nearly constant, and shows ~ 2 — 10 %. Since the SFR 
tightly correlates with LLya, and it roughly increases 
with the galaxy mass, some massive galaxies can be LAEs 
with Li^ya ^ 10^^ erg s~^. In addition, the /lae of LAEs 
having intrinsic Li^y^ ^ 10^^ erg s~^ change with cos- 
mic star formation history. However, the /esc decreases 
around the phase of SFR peak, and therefore suppresses 
the /lae. 

On the other hand, at lower redshift, the observa- 
tions indicate that n umber density of L AEs decreases 
by some factors (e.g., iCowie et al.ll201Ql ). The discrep- 
ancy may come from the difference in density field and 
the small box in our simulation. Our initial condi- 
tion is a somewhat special one which is focused on a 
MW-size galaxy, and the zoom-in simulation region is 
^ 5^ h~^Mpc^. Therefore, our simulation cannot re- 
produce the global statistics in observations. More- 
over, the LAEs fraction having EW > 25 A in this 
work shows ^ 44 % at z=4 and ^ 100 % at z = 6, 
which is somewhat h i gher t han th e LAE fraction in LEG 
sa mple (jStark et al.M2010lJ 2011: P entericci eFall 120111 : 
iSchenkeret alJ l l2012l : lOno et al . 20121 However, in obser- 
vation, the LAE fraction increases with decreasing UV 
brightness. Most of our model galaxies at z > 3 are 
fainter than the detection threshold in the LEG obser- 
vation. Since the number of galaxies brighter than the 
threshold of LEG observation is quite small (less than 
ten) , we need a larger sample covering a wide mass range 
to verify the model of LAEs. In addition, although some 
LAEs have been observed with UV continuum, and hence 
categorized as LEGs, it is inadequate to study LAEs from 
LEG-only sample, because a large fraction of LAEs may 
have UV continuum under the detection limit of current 
observations. We will address the general properties such 
as luminosity function, EW distribution and clustering 
systematically by using a set of uniform simulations with 
mean density field in larger volumes in future work. 
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5.3. The Viewing-angle Scatter of Escaping Lya 
Photons 

Despite their high metahicity, a fraction of galaxies at 
low redshift z ^ 1 show high escape fraction /esc of Lya 
photons (Figure [9]). We find that the escaping angle of 
the Lya photons depends strongly on the galaxy mor- 
phology and orientation, a phenomenon we dub as the 
"viewing-angle scatter" . Disky objects seen edge on can 
be hundred times fainter than the same objects seen face 
on. In a galaxy which has a gas disk, the Lya photons 
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Fig. 11. — The "viewing-angle scatter" of escaping Lyo; photons 
depends strongly on galaxy morphology and orientation. Shown 
here is the Lya escape probability in a irregular galaxy progenitor 
at z=3.1 (top panel) and the spiral MW galaxy at z=0 (bottom 
panel). The color bar indicates the probability per unit d(f),dcosO. 
There is no clear direction in the irregular galaxy, but in the disk 
galaxy, the Lyo; photons escape in a preferred direction normal to 
the disk. 



escape in a preferred direction normal to the disk, but 
there is no clear escaping direction in compact or irreg- 
ular galaxies without a gas disk. We demonstrate this 
effect in Figure [TTl We first estimate the normal direc- 
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Fig. 12. — The fraction of excitation cooling Lya to the total 
intrinsic Lya luminosity as a function of redshift. The red filled cir- 
cles indicate the values of the main progenitor, while the blue filled 
circles represent the median value of the galaxy sample indicated 
with black open circles. The error bars show the quart iles. 

tion to the gas disk according to the total angular mo- 
mentum of the gas, and set ^ = 0° along this direction. 
In a galaxy with irregular morphology such as the main 
progenitor at z = 3.1, there is no clear preferred escaping 
angle, as illustrated in the top panel of Figure [Til How- 
ever, in a spiral galaxy with rotationally supported gas 
disk such as the MW galaxy at z=0 in our simulation, the 
escaping angle is strongly confined to cosO :^ ±1, corre- 
sponding to 6> ~ 0° or 180°, as shown in the bottom panel 
of Figure [TTJ This is due to the fact that the Lya pho- 
tons have the minimum optical depth along the normal 
direction to the gas disk. More than 60 % Lya photons 
escapes to the direction of |cos 0\ < 0.5. Generally the 
Lya fiux from our model galaxies can scatter around the 
mean value typically by a factor of ten just from different 
orientations. 

As illustrated in Figure [3l most galaxies in our simula- 
tion have highly irregular shapes at high redshift due to 
accretion and gravitational interaction. At z=0, a num- 
ber of them evolve into spiral disks. The "viewing- angle 
scatter" explains why we see high Lya escape fractions 
in a number of low-z galaxies, and the fact that Lya is 
detected in a large numb er of face-on spira l galaxies in 
the nearby universe (e.g. JCowie et alll2010f ). 

5.4. Contribution of Excitation Lya Cooling 

There are two major mechanisms to generate Lya 
emission, the recombination of ionizing photons and the 
collisional excitation of hydrogen gas. However, the rel- 
ative contribution between the two mechanisms is not 
well understood. From our calculations, we find that the 
contributing fraction of excitation Lya emission to the 
total intrinsic Lya luminosity increases with redshift, as 
shown in figure [T2l 

In our cosmological simulation, galaxy evolution is ac- 
companied by cold, filamentary gas streams with tem- 



perature T ~ 10^~^ K, which penetrate deep inside 
dark matter halos (Zhu et al. in preparation, Ya- 
jima et al. in preparat ion), a ph enomenon also re- 
ported by other pcroups (jKatz et a l. 2003; Keres et al. 
2005, 2009; Birnboim & Dekel 2003; Dekel & Birnboim 
2006 : Q cvirk et al. 2008; Bro oks et al . 2009; Dekel et al. 
20091 ). Such cold gas can effic iently produce the ex- 
citation Lya cooling photons (iDiikstra fc Loebl l2009l : 
iFaucher-Giguere et al.l 120091 : iGoerdt et all I2010D . At 
higher redshift, galaxies form through more efficient gas 
accretion and more frequent merging event. As a re- 
sult, the contributing fraction increases with redshift, 
and becomes dominant at z > 6. This excitation mech- 
anism does not depend on the stellar radiation, and can 
therefore produce high Lya EWs. We find that the 
EWs of LAEs increases significantly at z > 6, reaching 
>10^Aatz~10. This is larger than the upper- limit of 
EW, 240 A, which considers only stellar sources assum- 
ing a Salpeter IMF with solar abundance of metallic- 
ity ([Chariot fc Falllll993D . Although the upper-limit in- 
creases with decreasing metallicity, it was suggested that 
top-heavy IMF like Pop III stars are needed for making 
EW > 400 A (|Schaererl 120031 : iRaiter et al]l2010D . How- 
ever, even though Salpeter-IMF is used in this work and 
the stellar metallicity is mostly Z/Zq > 10~^, we find 
that the EW can be higher than the upper-limit by the 
efficient excitation Lya emission. On the other hand, the 
Lya line i s stro ngly damped by IGM correction at z > 6 
(iHaimanl 120021 : iLaursen et 311120111 : iDaval fc Libeskindl 
I2OIII ). which can result in a lower EW. The suppression 
by IGM highly depends on the i nhomogeneous ioniza- 
tion structure around LAEs (e.g., [M cQuinn et al.l 120071 : 
Mesinger fc Furlanettol 120081 : llliev et al. 20Q8[) . We will 
address the detectability of high-redshift LAEs and EW 
after IGM correction by running large-scale Lya RT in 
IGM in future work. 

5.5. Lya Luminosity Functions 

The simulation box in this work is too small to study 
global statistics directly. As a rough estimate, we use 
the luminosity - halo mass correlation we find above may 
be used to construct Lya luminosity functions (LFs) at 
different redshift when combined with halo mass func- 
tions from large-box, general cosmological simulations. 
For example, at z = 3.1, we divide all galaxies in the 
snap shot by the halo mass with 0.25 dex, and fit to 
the median value of each bin, this gives a correlation of 
LLya (erg s-1) = 10^^-^^ X (MtJ^/Mc^) We then use this 
to con vert the halo mass function of iSheth fc Tormen I 
(|1999[ ) to the Lya LF. 

Figure [13] shows the resulting Lya LFs in comparison 
with observations at redshift z = 3.1,5.7, respectively. 
The red solid curves are LFs above a detecti on threshold 
of S-[ ^ ^fy = 10~^^ ergs s~^ cm~^ arcsec"^ (jOuchi et al.l 
I2OO8I . I2OIOI ). while the red dashed lines represent LFs 
from total luminosity (counting all escaped photons with- 
out a fiux cut). 

While the un-filtered LFs se em to agree with observa - 
tions of lGronwah et all (|2007D and lOuchi et all (|2008r ). 
the filtered ones are significantly off. The difference 
comes from the reduction of Li^y^ and /esc due to the 
fiux cut. Moreover, the dispersion in the luminosity - 
halo mass relation at different redshift may cause a large 
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Fig. 13. — Derived hya luminosity functions from our 
model in comparison with observations. The red, solid lines 
are filtered LFs above a detection threshold of S'Lya = 
lQ-18 gj.gg g-i cm~2 arcsec"^ (jOuchi et al.l l2008'). while the red, 
dashed lines represent unfiltered ones. The black dashed lines are 
the Schechter function at z = 3.1 and 5.7 derived from LAEs ob- 
servation by Ouchi et al. (200 ^ ), and the black dotted line is from 
observation at z = 3.1 bv lGronwall et aTl (|2007l ). 

scatter in the LFs. This plot suggests that the current 
simulation in this work is not suitable to study a large 
galaxy population and its statistical properties, because 
there are too few observable LAEs. Moreover, as dis- 
cussed earlier, the predicted Lya properties may be af- 
fected by a number of numerical and physical limitations 
of our model. For example, the one-phase model cur- 
rently used in the present work may underestimates the 
density of cold hydrogen gas, and hence underestimates 
the Lya flux. We will study the Lya LFs at different 
redshift in a forthcoming paper with the improved ART^ 
which incorporates a two-phase ISM model, and a gen- 
eral simulation with mean overdensity in a larger volume 
(Yajima et al, in preparation). 

5.6. Limitations of Our Model 

As demonstrated above, our model is able to explain 
a number of observed properties of LAEs at different 
redshift. However, we should point out that our current 
simulations suffer from a number of major limitations 
which may affect the predicted Lya properties. 

• In the current work, we use a one-phase ISM model, 
which considers the average density and tempera- 
ture of the gas. Such a model likely underestimates 
the density of cold hydrogen gas, which may lead 
to significant underestimate of the Lya emission 
coming from cold 10^ K), dense gas. On the 
other hand, such a model also underestimates the 
amount of dust associated with cold molecular gas, 
which likely results in underestimate of absorption 
of Lya photons by gas and dust. We will inves- 
tigate the Lya RT and ionization structures in a 
two-phase ISM model in a forthcoming paper. 



• The absorption and transmission of IGM are not 
taken into account in the RT calculations. As dis- 
cussed in the previous section, these two effects 
may suppress the Lya flux and change the line pro- 
files. 

• The simulations do not have sufficient resolutions 
to resolve dense regions and outflow, w hich requires 
a hig h spatial resolution of ~ pc (e.g.. iFujita et al.l 
[20091 ). It is a challenge for cosmological simulations 
to resolve both the inflow gas from large scales of 
~ Mpc and the outflow from pc-scale star forming 
regions. For a simulation with a box of 100 Mpc 
like the one we have, this requires a large dynami- 
cal range over eight orders of magnitude, which is 
beyond the scope of our current work. 

• The simulation box in this work is too small to 
study a large galaxy population, as well as effects 
of environment on galaxy properties and their evo- 
lution. One needs uniform simulations in large vol- 
umes in order to systematically investigate the for- 
mation and evolution of L* galaxies. 

Finally, we stress once again that caution should be 
taken when comparing directly the results from our cal- 
culations to data from a given survey, because, as dis- 
cussed above, the observed Lya properties depend sen- 
sitively on a number of factors, including galaxy proper- 
ties, viewing angle, model parameters, and observational 
threshold. 

6. SUMMARY 

To summarize, we have investigated the Lya properties 
of progenitors of a local L* galaxy by combining cosmo- 
logical hydrodynamic simulations with three-dimensional 
radiative transfer calculations usine: the new ART^ code. 
Our cosmological simulation follows the formation and 
evolution of a Milky Way-size galaxy and its substruc- 
tures from redshift z — 127 to z = 0. It includes im- 
portant physics of dark matter, gas dynamics, star for- 
mation, black hole growth, and feedback processes, and 
has high spatial and mass resolutions to resolve a MW- 
like galaxy at z=0 and its progenitors at higher redshifts. 
Our radiative transfer couples Lya line, ionization of neu- 
tral hydrogen, and mult i- wavelength continuum radiative 
transfer, which enables a self-consistent and accurate cal- 
culation of the Lya properties in galaxies. 

We find that the main progenitor of the MW galaxy 
is Lya bright at high redshift, with the emergent Lya 
luminosity close to the observed characteristic ^Lya 
LAEs at z ~ 2 — 6. Most of the fainter galaxies in the 
simulation fall below the detection threshold of many 
current surveys. The Lya escape fraction correlates with 
a number of physical properties of the galaxy, such as 
mass, SFR and metallicity. We find a "viewing- angle 
scatter" in which the photon escape depends strongly on 
the galaxy morphology and orientation, such that the 
Lya photons escape in a preferred direction normal to 
the gas disk in disk galaxies, but randomly in compact or 
irregular galaxies. Moreover, the EWs of LAEs increases 
with redshift, from tens of Angstroms at redshift z ~ 
to hundreds of Angstroms at z ~ 8.5. Furthermore, 
we find that high-redshift LAEs show Lya line profiles 
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characteristic of gas inflow, and that the Lya emission 
by excitation coohng increases with redshift, accounting 
- 50 % of the total at z > 6. 

Our results suggest that galaxies at high redshift form 
through accretion of cold gas, which accounts for the the 
high EWs, the blue-shifted line profiles, and the domi- 
nant contribution from excitation cooling in Lya emis- 
sion. Moreover, some of the observed LAEs Sit z ^ 2 — 6 
with I/Lya ^ io^^~^^ erg s~^ may evolve into present-day 
L* galaxies such as the Milky Way. 

We thank Carlos Frenk for kindly providing us the 
Aquila initial condition for the cosmological simula- 
tion. We thank Mark Dijkstra, Claude- Andre Faucher- 
Giguere, Eric Gawiser, Lars Hernquist, and Avi Loeb for 
stimulating discussions and helpful comments, as well as 
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